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1. Introduction 

Many interesting observables in hadron physics, e.g. (moments of) generalised parton distri- 
butions (GPDs) or distribution amplitudes (DAs), can be computed from matrix elements of local 
operators between hadron states. (For an extensive review of GPDs see Ref. [Hp, for DAs see, 
e.g., Ref. [Q].) GPDs parametrise a large class of hadronic correlators, including form factors and 
the ordinary parton distribution functions. Thus those distributions provide a formal framework to 
connect information from various inclusive, semi-inclusive and exclusive reactions. Furthermore 
they give access to physical quantities which cannot be directly determined in experiments, like 
e.g. the orbital angular momentum of quarks and gluons in a nucleon (in a given scheme) and the 
spatial distribution of the energy or spin density of a fast moving hadron in the transverse plane. 
Since the structure of GPDs is rather complicated direct experimental access is limited. Therefore, 
complementary information from lattice QCD should be used. For a recent overview on Monte 
Carlo results, see Ref. [^J. 

To make contact with the continuum, the lattice operators needed for GPDs and DAs have to 
be renormalised. Compared with moments of ordinary parton distributions, the specific difficulty 
in the treatment of moments of GPDs and DAs lies in the fact that the required matrix elements 
are no longer forward matrix elements. This circumstance complicates the pattern of mixing under 
renormalisation. 

In Ref. [|J] we have calculated the non-forward quark matrix elements needed for the renor- 
malisation of quark-antiquark operators with two derivatives, which determine the second moments 
of GPDs, and we have discussed the mixing problem in detail. This calculation was performed in 
one-loop lattice perturbation theory for the Wilson fermion action. 

Here we present some of our results [||] obtained for 0(a) improved fermions using the 
Sheikholeslami-Wohlert (clover) action ^ (without operator improvement) and new calculations 
for overlap fermions p]. 

The fermion action is of the following generic form (spinor and colour indices are suppressed) 

S F = yD\ l / = a 4 Y,9(x)D(x-yMy) (1.1) 

x,y 

with 

D sw = D w - i a g r c sw - [y„ , y v ] F^° veT , (1.2) 

D N = =P-(l + -?L=) , X = D W ~- (1.3) 

for the clover (Dsw) an d overlap (Djv) operator with massless quarks, respectively. Here Dw de- 
notes the Wilson-Dirac fermion operator with the forward (backward) covariant derivative A^ (At) 

D w = \ [YLi( A li+ A n)- arA li A n] i a mV(*) = - [U Xi n\i/(x+afl)-\i/(x)) . (1.4) 

g is the bare coupling, a the lattice spacing and F^° veT the standard "clover-leaf" form of the lattice 
field strength 

f dover = 1 y (tjD _ ^Df \ q 5) 
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In the perturbative calculation the operators to be investigated are sandwiched between off- 
shell quark states with 4-momenta p and p'. Our calculations are performed in a general covariant 
gauge, the final numbers are presented in Feynman gauge for the Wilson parameter r = 1 and the 
overlap parameter p = 1.4, leaving the coupling strength of the improvement term c sw free. 

In [Q] a detailed discussion of the renormalisation procedure for the case of non-vanishing 
momentum transfer can be found. The matrix of renormalisation and mixing coefficients Zn(a\l) 
relating regularised lattice vertex functions T^(p' ,p,a,gR) and MS renormalised vertex functions 
Tf{p\p,jJL,g R ) is defined such that 

N 

rf(p',p,n,g R )=Z ¥ l (an) J £z ik (an)r L k (p\p,a,g R ) (1.6) 

k=\ 

with the quark wave function renormalisation constant Z v . Here p (//) denotes the momentum of 
the incoming (outgoing) quark, the renormalisation scale is fi, the renormalised coupling is denoted 
by gR, and N is the number of operators which mix in the one-loop approximation. 

2. Operators and mixing 

We consider quark-antiquark operators with up to two covariant symmetric lattice derivatives 
D = D — D (for their definitions, see e.g. []|]) and external ordinary derivatives d. 

In the non-forward case q = p' — p ^ we use two realisations of operators with covariant 
derivatives in momentum space: either the momentum transfer "acts" at the position x associated 
with the operator (realisation I) or q is applied at the point half way between the quark fields iff and 
Y (II)- As an example, for two covariant derivatives this leads to 

«- « \(«) 1 „ ._ 

YD^D v y) {q) = —2^{¥(x)U x ^U x+afl , v \l/(x + afl + av) 



-y(x + av)U x+a $ >tl ul +afL ^(x + afL) 
- \jr(x + afi)u£„ U xy y(x + av) 

+ W + ^+aV)Ul +a , il UlMx)) e ^+afl/2 + av/2) ^ 



and realisation I is obtained from (2.1) as 



yD^Dvyj W (q) = cos (^) cos (^) (yD^D^ W (q) . (2.2) 

We have set the Dirac matrix in the operators equal to the unit matrix for simplicity. 

In what follows, we indicate in the operator symbols the derivatives by superscripts D and d. 
The quark-antiquark operators with one derivative are given by 

<v = - \ YYnDy ¥ , <v° = ~9Yn Y&v ¥ , (2-3) 



f lvco = -~¥[Yii,rv]Da¥, K'vco = ~~dw {W^YvW) ■ (2-4) 
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T D 

Operators such as ^a'va are of interest for tensor GPDs as well as for transversity and we call them 



transversity operators. For non-chiral fermions, operators ( |2.4| ) contribute as lower-dimensional 
operators to mixing in certain representations which determine second moments of GPDs. 
As operators with two derivatives we consider 



1 



• M vo>- jVYiiDyDcoY, @lva = -j d AvYiiD< w)i 0™ a = -jd v d a [WYtiW) , (2-5) 
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— c 
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1 



»5,DD 
•71VCO 



^YYtiYsDvDay, = ~ jd v (Vy^D^) , 6>$l = --d v d a {WnYsV) 



and the transversity operators 



^vwa = -^[YiiilfrPaDoY, tfivia = -~d m d a {W^Yv\w) 



1 

V 



(2.6) 



On the lattice the operators are classified according to the irreducible representations Zjr of 
the hypercubic group H(4) (for the notation see, e.g., Ref. Here / denotes the dimension of 
the representation and k labels inequivalent representations of the same dimension. In addition our 
operators will be chosen such that they have definite charge conjugation parity C. 

To define various index combinations of operators we use the shorthand notations: 



^-{v 1 v 2 } 

^{VlV 2 V 3 } 

^ > ||v 1 v 2 v 3 || 
'((V1V2V3}) 



•V1V2 



+ - 



■V 2 V l J 1 



" l^Vl V 2 V 3 4" V ViV3V 2 4" ^V 2 Vl V3 4 
^ViV 2 V3 — ^Vl V3V2 4 ^V3ViV 2 — ^V3V 2 V[ 
^VlV 2 V3 4 ^ViV3V 2 — &Vt,V[V2 ~ &V}V2Vl 



V 2 V3V[ 4 

21 



&V3V1V2 4 ^V3V 2 Vi ) ) 
? V 2 V3Vi 42^ v 2 V[ V3 J 



(2.7) 



For the first moments we choose the following representations and operators: 



Operator 



fD 
{14} 



,5,D 
{14} 

5.D _ 1 
44 3 
T,D 

{{124» 
T,D 

((122)) 



Representation C 



? T,D 
((133)) 



? 5,D 
33 



,(«) 
'3 

r (3) 
'1 

(6) 
'4 

r (3) 
'4 

'2 

,(8) 



(2.8) 



All operators in (2.8) are multiplicatively renormalisable. These representations exhaust all possi- 
bilities for the twist-2 sector in the continuum. 
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In these proceedings we present results for the following sets of mixing twist-2 operators 
which are related to the second moments of GPDs in the unpolarised case: 



c 



1 



r (8) 



c 



5>DD 
{114}" 

((114)) ■ 

? 5,dD 

1121311 ' 



'{124} ' 



sDD 



5>33 
{124} 



{224} 



+ 



{334} ) 



{"4} 



'{224} 



'((224)) 



+ 



7 ((334)) 



f 5.dD 

{{213)) ' 



f 5,DD 
12131 



In the polarised case we consider 



j 4) ,C = +l : 
T ? (8) ,C = +1: 



? 5,DD 
{124} ' 



{124} 





\ 


- ^{114} 




/?5,DD 


1 


- ^«114» 


~ 2 


- „ 


G\ 


213 ' 





+ 



f 5,DD 
{224} 

V ((ZU}) 



+ 



{{213)) ' 



{{114)) " 

- ^ r ' 5 - 

- ^ 411 



1 

2 
1 

" 2 



? 33 
{{224)) 

y 422 



5,33 \ 

{334} ) 



+ 



{{334)) 



+ 



? r,3 

433 



{334} J ' R 
{{334)) ) 



{114} 



f DD 

1 1 2 1 3 1 1 ' 



1 

2 

5,33 
{{114))" 

^123 



{224} 



+ 



=,5,33 \ 
{334} J 



z®5,33 
^((224)) 

? 23? _ ' 



+ 



5,5,33 \ 

{{334))J 



5,T,D 
132 • 



(2-9) 



(2-10) 



(2-11) 



(2-12) 



3. Results 

We calculate the matrix elements in one-loop lattice perturbation theory in the infinite volume 



limit following Kawai et al. [10]. The renormalisation factors for operators related to the first 
moments of GPDs are the same as in the forward case and known for the popular representations. 

In the case of the second moments of GPDs and distribution amplitudes we present the renor- 
malisation matrices in the form 



Z ik (an) = 8 ik -^^(y ik \n{a 2 pL 2 )+B ik ) 



(3.1) 



where 7 denotes the matrix of anomalous dimensions and B is the matrix of finite one-loop contri- 
butions. Below, the numbers in boldface indicate results known from the forward case. 
In the representation < $L% two operators mix in the non-forward case leading to 




B 



overlap 



-47.4441 -0.95719 \ 
-17.418 J 



(3-2) 



B 



clover 



-11.563 + 2.898c™, -0.9844, 0.024 -0.255c sw -0.016c? 







20.618 + 4.746c w -0.543ci 



(3-3) 



The results do not depend on the particular choice of the covariant derivative in the non-forward 
case. 
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For the representation (2.10) six operators of the same dimension mix (in the one-loop ap- 
proximation the contribution of G7 is absent). The finite contributions for clover fermions B clover = 

+ c sw B^ + c 2 w B^ can be found in ^ for the two different derivatives. For the choice (I) of 
the derivative we find for overlap fermions 
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5 
6 
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7 
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5 
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2 
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-2 


\o 
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I / 



(3.4) 



^overlap 



V 



48.1089 



3.31602 







-3.43393 
-17.4180 
20.2671 


-12.2181 
-12.2181 



0.53795 



46.8416 







1.33526 


-10.9597 
17.41796 





0.03132 


-3.52389 


-34.1678 
0.6396 



0.64459 \ 


1.79504 


1.91865 
-32.8888 / 



(3.5) 



Using clover fermions and the representation (2.10), there is a dangerous mixing between the 
operator ff\ and the operator which is one dimension lower than G\\ 



^clover I 

1 \\/a— part 



8 2 rCf 
I6n 2 



(-0.51771 +0.08325 c JVV -0.00983 



1 



sBorn 



(3.6) 



This mixing leads to a contribution which diverges like the inverse lattice spacing in the continuum 
limit. Thus the perturbative calculation of the mixing coefficient is not reliable and the operator 6% 
has to be subtracted non-perturbatively from G\. Note, however, that G\ and G% are of opposite 
chimlity. 

Using chiral lattice fermions in perturbation theory, the general matrix element of an operator 
with two co variant derivatives contains also 1 /a contributions indicating the potential presence 
of lower dimensional operators at one-loop level. Choosing a particular irreducible representation 
such as ( 2.10| ), the contributing terms are expected to cancel. This indeed happens in our calculation 
with overlap fermions and the mixing is absent to a very high accuracy numerically: 



^overlap i 



\ja— part 



0. 



(3.7) 



Using the operators %L. 1 1| ) and ( J2. 12[ ) (in the polarised case) the one-loop results for the chiral 
overlap fermions coincide with those given in (EO) and (Oh,(|0), respectively. 



4. Summary 

We have considered quark-antiquark operators needed for the computation of the first two 
moments of GPDs and meson distribution amplitudes within the framework of lattice QCD. In 
one-loop lattice perturbation theory we have calculated the non-forward quark matrix elements of 
these operators employing clover improved Wilson fermions and overlap fermions together with 
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Wilson's plaquette action for the gauge fields. From the results we have determined the matrices 
of renormalisation and mixing coefficients in the MS-scheme for some sets of operators belonging 
to particular representations of the hypercubic group. 

If there is only mixing between one operator with two covariant derivatives D and one oper- 
ator with two external derivatives d (e.g. for the representation (23) with three different Lorentz 
indices) the mixing coefficient turns out to be rather small. In the case of the representation ( 2.1C ) 
with two equal Lorentz indices - which contains eight potentially mixing operators - the mixing 
effects are more severe. As expected, the dangerous mixing with a lower-dimensional operator of 
opposite chirality in the case of clover or Wilson fermions is absent using overlap fermions. This 
obviously favours the use of chiral fermions in future simulations related to GPDs and DAs. 

Results for a more complete set of representations using overlap fermions will be presented 
elsewhere [11]. 
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